Boundary induced amplification and nonlinear instability of interchange modes 
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It is shown that small distortions on the boundaries are amplified in the core of a magnetized 
plasma if the system is close to marginal stability for the ideal magnetohydrodynamic interchange 
mode. It is also shown that such marginal systems can be nonlinearly unstable. The combination of 
boundary amplification and nonlinearity is shown to result in a nonlinear instability. The induced 
instability is highly sensitive to the boundary in that, if the fractional deviation from marginality 
is a small parameter b, the system can go unstable from fractional boundary distortions of 0(b 3 ^ 2 ). 



I. Introduction. — Magnetically confined plasmas for fu- 
sion are limited in how much pressure can be contained, 
by the so-called f3 limit - the critical ratio of pressure 
to magnetic energy density. [T] This limitation generally 
comes from interchange instabilities, wherein flux tubes 
of high pressure plasma can interchange with outer, lower 
pressure flux tubes. [XH3] Such energy release can be stabi- 
lized if the flux tube interchange is disallowed by topology 
(on account of the frozen-in condition for strongly mag- 
netized plasmas), that is to say, if the "transverse" field 
is strong enough. [3] For maximum efficiency, one wants 
to operate close to marginal stability j3 — > j3 c . 

In this letter, we report on two findings pertaining to 
operating near marginal stability. Let A/3 — /3 C — /3 and 
let the system size be a. Then, (1) we establish the gen- 
eral idea that a small perturbation of S /a on the bound- 
ary amplifies interchange displacements in the core of 
the plasma, by an amplification factor /3 c /A/3; (2) we 
show that the system is nonlinearly unstable, so that a 
critical boundary perturbation will destabilize the inter- 
change mode even for systems below the linear /3 stability 
limit. Upon combining these two findings, we find that 
the amplification phenomenon leads to a nonlinear in- 
stability criterion which is highly sensitive to boundary 
perturbations, namely the fractional critical size is even 
smaller than the fractional deviation from marginality, 
i.e., (S/a) > | A/3//3 c | 3 / 2 . This has the implication that 
magnetic configurations designed to confine plasma close 
to the j3 limit within a tolerance of e would necessitate 
that the design be more sensitive to boundary perturba- 
tions; specifically, boundary tolerances need to be better 
than e 3 / 2 . Such considerations are of significant impor- 
tance in the design of axisymmetric tolerances for ad- 
vanced tokamaks as well as in the fully 3D design of stel- 
larators for fusion. 

Field amplification near marginal stability has been 
shown to also occur in tokamak plasmas, for kinklikc[4 
and tearing[Sl |B] modes. It was reported in Ref. [5] that 
external perturbations would be amplified if the equilib- 
rium profile is close to marginal stability for tearing. In 
Refs.|Hand|6l it was shown that error fields at the plasma 
edge are amplified in the core by a factor that is inversely 
proportional to the marginal stability parameter. Our 



FIG. 1. The conducting plate boundaries, at x = ±a, of a 
system with Vp = p x opposite a gravitational force g = — gx 
balanced by a transverse field Bo = Soy are perturbed by a 
ripple of amplitude 8. 

present result shows that the amplification phenomenon 
extends also to interchange modes (as also found in a 
related case in Ref. [7]). We include nonlinear perturba- 
tions in our methodology to show that the amplification 
precipitates an already dormant nonlinear instability. 

We use a model system to illustrate the basic phe- 
nomena. The calculation is done in slab geometry with 
an effective gravitational field to model field line curva- 
ture. Boundary perturbations are treated as a ripple on 
the boundary conditions, similar to the Kulsrud-Hahm 
problem. [5j 

II. Linear boundary perturbation problem. — We be- 
gin with a simple-minded linear calculation to first 
demonstrate the amplification phenomenon. Con- 
sider an incompressible system describable by the 
two-dimensional, reduced magnetohydrodynamic (MHD) 
equations [3] given by, 

dt1> = {il>M> (!) 
z.ViX^u + ^u}) = W,V^} + S 9 s p, (2) 

u = zxV|_<^, Bj_ = zx VlVj (3) 
where, in general, p = p(ip) and 

{f,h} = d X fdyh-d y fd X h. (4) 
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The problem of interest is illustrated in Fig. [T] with 
p' and Bq constant, and such that we are stable to the 
ideal MHD interchange mode. The system is taken to be 
periodic in the y direction, and 



d y ip = —ksin(ky)Sd x ip 



(5) 



at x = ±a — Scoa(ky) provides the rippled boundary 
condition in the x direction. For 5 = this system is in 
static equilibrium. 

We now perturb the boundaries at x = ±a by adiabat- 
ically introducing S/a <C 1 and allowing for a quasistatic 
equilibrium, ip = Box + ip, to form. In general, ip satisfies 
the nonlinear equation 



and the boundary condition, Eq. (|5j). The solution 
6B 



ip = — " ; cos(k x x) cos(ky), 
cos(k x a) 



where 



, 2 _ 9P0 _ , 2 



(6) 



(7) 



(8) 



satisfies Eq. (6| and Eq. ([5| to lowest, linear, order in S/a. 
From Eq. (17), we note that as k x a approaches tt/2, the 
perturbation gets amplified. In fact, as we will confirm 
later, k x a = tt/2 corresponds precisely to the linear sta- 
bility criterion for the ideal interchange mode, where the 
critical transverse field for marginal stability, Bq = B c , 
is obtained from Eq. ^ at k x = k c = it /2a. 

We can determine the scaling of this amplification by 
letting Bq = B c + b and k x = k c — Ak xi where b/B c ~ 
Ak x /k c <C 1. To lowest order in b/B c 



Ak = 



where k\ = k 2 + k 2 , and Eq. (|7| yields 



S/a k c 
b/B< 



k\ 



B c cos(k x x) cos(fcy). 



(9) 



(10) 



Therefore, as we approach criticality by letting b — > 0, a 
small perturbation at the boundary of 0(S/a) can induce 
a large response in the bulk of the plasma, scaling like 
B c S/ba. At exactly Bq = B c the solution becomes ill- 
defined away from the boundary. 

The solution given by Eq. ^ is valid as long as S/a 
is small, such that B C S ba. However, this solution 
clearly establishes the phenomenon of amplification. We 
show in what follows that nonlinear effects arise at even 
smaller boundary amplitudes, significantly modifying our 
understanding of the effect of boundary perturbations on 
stability. 



777. Nonlinear evolution. — For a system that is 
marginally stable to the ideal MHD interchange mode, 
it was previously shown that small homogeneous pertur- 
bations in the plasma can result in nonlinear, explosive 
growth. [TU] The nonlinear instability result showed that 
the optimal scaling of magnetic perturbations was given 
as IV'/V'ol ~ e = {b/Bc) 1 / 2 . This realization prompts 
us to apply this nonlinear stability scaling as optimal 
ordering for nonlinearities in the present amplification 
calculation. Thus, optimally, we should have S/a ~ e 3 , 



using Eq. (10). With this scaling, we will show that S 



on the boundary can introduce a nonlinear instability in 
the plasma. It is important to note that even though we 
order the parameters as described, the boundary pertur- 
bation amplitude, S/a, and marginality condition, b/B c , 
are independent, small parameters. 

Using the marginality condition as a smallness param- 
eter we expand ip in a series, i.e. let 



(11) 



where each successive term is smaller by a factor of e and 
ipo = (B c + b)x. By matching terms order by order, we 
solve Eq. |6]), where ^ = ip — ipQ, using the boundary 
condition given by Eq. 
To order e, Eq. ^ yields 



B 2 c {Wl + fclKVi = 0, 



(12) 



where we have substituted for gp' using Eq. Q with 
Bq = B c and k x = k c . Taking Eq. ^ to lowest order 
implies that 



= ±a 



= 0. 



(13) 



Including the lowest order term in the solution given by 
Eq. 0, we find that 



A 



S/a k c 



B, 



i(k c x)cos(ky), (14) 



where, for convenience, A is introduced as a free param- 
eter in this particular manner to represent the plasma 
response. We are interested in how S in the boundary 
induces A in the plasma. We use the scaling k c A/B c ~ e 
for the plasma perturbation, but it will be treated as a 
separate small parameter. 



The lowest order equation, given by Eq. (12), can also 



be arrived at by solving the linear, ideal MHD, inter- 
change mode problem and insisting that lo = 0. This 



results in a zero frequency state, given by Eq. ( 14 ) with 



(5 = 0, where Alfvenic restoring forces exactly balances 
the Rayleigh- Taylor growth rate, i.e. k\y\ c — gp' /po = 0. 
From this, we can conclude that B c is the critical field 
strength needed to be at marginal conditions. Allowing 
for b > means that the system is marginally stable to 
the ideal MHD interchange mode. 
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Continuing to order e 2 , we find that ip2 satisfies the 
equation 



B 2 (V 2 +kl)d v ^ 2 



-£ c {V>i,V 2 tM (15) 
= (16) 



and the boundary condition 



0. 



(17) 



This implies that ?p2 is only a function of x. We can solve 
for -02 by taking Eq. ([!]) to second order and averaging 
over y to get 



(18) 



where the bar denotes an average over y. Equation ([!]) 
to first order implies that 



d t ipi = B c d y (p, 



and so we find that 



4'2 



1 fee 

4: B r 



A 



S/a k c 
b/B c k? ± 



B, 



(19) 



sin(2fc c x). (20) 



This term represents the flattening of the perturbed field 
lines (zonal field) to second order, driven by the first 
order perturbation. 

In order to find how 5 drives the amplitude A we extend 
our analysis to order e 3 where Eq. |6| yields 

B 2 C {V\ + k 2 ± )d y yj 3 + 2B c b^\d y ^ = 

-B e ({^,W 2 M + {iJ2,V 2 M), (21) 

with the boundary condition 

d y ^\ x=±a = -SB c ksm(ky). (22) 

The above equations imply that ip3 = ip3(x) cos(ky) 
where ip3(x) is a linear combination of k c and 3fc c har- 
monic terms. The 3fc c harmonic is straightforward and 
uninteresting; however, the k c harmonic is secular, so we 
will focus on obviating this secularity. Using Eq. Q to 
order e 3 we let 



1p3 



B c —x sin(fc c x) cos(fcy). 
a 



(23) 



Substituting Eqs. (14), (20), and (23) into Eq. (21) and 



insisting that the secular terms go to zero yields 
k 2 k 2 - 3k 2 r , S/a k 



2B c bA 



.4 



b/B c k 



B r 



0, (24) 



after simplification. The above result gives the sought- 
after relationship between the boundary perturbation 
amplitude and the amplitude of the plasma response. 

To discuss this result we now consider allowing for the 
plasma response to evolve in time, at a rate slower than 
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FIG. 2. Plot of the potential energy U(A; S) as a function of 
amplitude, A, for: (a) 5 = 0; (b) 5 < 5 C ; and (c) S > S c . The 
dotted box shows the shrinking boundaries of the stable well 
(shaded region). 



Alfven time, ta = {k c VA c )~ l , as we distort the boundary 
at an even slower rate. Explicitly, we let A = A(t) with 
T~Adt ~ e but we keep 5 small. It is easy to show that 
doing this results in the nonlinear time evolution of A(t) 
given by 



A = -2bA 



l fc 2 -3 
4/c 2 + l 



A 



2 <5 1 
7T b k 2 + 1 



(25) 



where we normalize the variables by setting k r . B r , and 
Vac equal to 1, for simplicity. Multiplying Eq. (25) by A 
and integrating once yields the "energy" integral E Q = 
A 2 /2k 2 + U(A;S), where 



U(A; 5) = bA 2 - 



1 k 2 



16 k 2 + 1 



A 



2 5 1 



I 4 



7T b k 2 + 1 



(26) 



The above equation represents a potential energy as a 
function of A, given k, b and S, and determines the overall 
stability of the system. We have chosen b > 0, to be 
marginally stable, and so, for k 2 < 3 the system is stable 
to all perturbations. However, for k 2 > 3 and fixed b the 
stability of the system is dependent on the size of S and 
A = A(t = 0). " 

The potential, with k 2 > 3 and fixed b, is shown in 
Fig. [2] for different values of i5. For 5 — (Fig. [2^,), there 
is a stable well for \A \ < A c , where 



A c = 2y/2 



k 2 + 1 



k 2 



b 1 ' 2 . 



(27) 



This result is the same as the one found in Ref. 10 where 
nonlinear stability was noted and simulated for large 
enough A. As 5 is increased the stable well shrinks as 
the two positive roots of U'(A;S) merge. From Fig. J2Jd 
we can see that the symmetry of the potential is broken 
and one side of the stable well drops so that the critical 
Aq to stay nonlinearly stable is less than the one given 
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in Eq. (p7|. When 6 = 5 C , where 



Sr.= 



the positive roots of U'(A; S) become degenerate and the 
stable well becomes a point. When 5 > 8 C (Fig. [2j;) the 
system is always unstable for any Aq. This means that, 
even though the system is linearly stable, a boundary per- 
turbation of order (b/B c ) 3 / 2 can precipitate a response in 
the plasma of much larger amplitude, of order (b/B c Y' 2 , 
accompanied by explosive growth. 

We remark that, in this calculation, we have assumed 
the density gradient, p' , to be constant. This is done for 
simplicity, to illustrate the two phenomena of amplifica- 
tion and nonlinear instability in a transparent manner. 
We have shown elsewhere that keeping a general po(x) 
did not result in new phenomena or mitigate the appear- 
ance of the foregoing phenomena, and that p' constant 
is a reasonable simplification. 

IV. Conclusion. — We have shown that for systems op- 
crating close to marginal stability for the interchange 
mode, a small perturbation on the boundary can induce a 
large response in the core of the plasma. A simple linear 
analysis shows that a secondary equilibrium, with ampli- 
tude inversely proportional to the marginal stability pa- 
rameter, can exist as boundary perturbations are added 
adiabatically. This suggests that the system would need 
to be far enough away from criticality so as to stay well- 
defined. We accordingly extend this analysis nonlinearly 
to show that, even if the perturbation were scaled much 
smaller than the marginal stability parameter, the sys- 



tem can become nonlinearly unstable for boundary per- 
turbations larger than a critical value. These results have 
implications in the stability analysis and corresponding 
design of magnetic confinement devices operating close 
to marginal conditions. With A/3 > 0, a linear stability 
analysis would show that the system is stable to all per- 
turbations; however, relatively much smaller boundary 
perturbations, of order (A/3//3 c ) 3 / 2 , can destabilize the 
system nonlinearly, with subsequent perturbations grow- 
ing without small amplitude saturation. 
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